Theory of Upper Critical Field without Energy Quantization 
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Conventional theories for determining upper critical fields are inevitably related to the lowest 
eigenvalues of appropriate equations. In this Letter, a new theory of upper critical fields is designed 
and justified. Using MgE>2 as modeling prototype, our computations are in excellent agreement with 
the Ginzburg-Landau theory. The long-standing issue, the upward curvature of the upper critical 
field, is found to be a manifestation of the crossover of the order parameter. The current theory 
is an alternative to the traditional technique of energy quantization in determining upper critical 
fields. 
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The upper critical field {B C 2) is the maximum magnetic 
field that a bulk type-II superconductor [lj can sustain 
in the superconducting state. The value of B C 2 is very 
important as it partially determines the current carrying 
capacity of the superconductor and its uses e.g. to pro- 
duce high field superconducting magnets. Furthermore, 
study of properties of B c2 may test the validity of various 
theoretical models and provide information for impor- 
tant superconducting parameters such as the coherence 
length. Hence, research of the upper critical field is of 
practical, fundamental and enduring interest. 

To determine the upper critical field, an important 
starting point is the linearized Ginzburg-Landau (GL) 
equation |2j and its variants. Examples would be the 
quantum harmonic oscillator equation hi, the Mathieu 
equation |3(, the p-wave GL equation and the d- 
wave GL equation eJ. Determining upper critical fields 
from these equations involves finding the lowest eigen- 
values of the said equations. Alternatively, the micro- 
scopic description of the upper critical field, based on 
the Gor'kov gap equation [6(, may also be reduced to 
findingthe lowest eigenvalues of appropriate equations 
0H B E3- Note that the eigenstate Q of the low- 
est eigenvalue of the linearized GL equation is usually 
chosen as the trial function in the variational deter- 
minations of B C 2 (see, e.g., Refs. 0, EI). Thus, the 
variational treatment is related to the lowest-eigenvalue 
method. Note further that the perturbation approach 
to determining upper critical fields is also intimately re- 
lated to the scenario of the lowest eigenvalue Conse- 
quently, we are led to conclude that almost all the efforts 
El II 11 II II II 0, IE i, ES IH E2 for determining up- 
per critical fields are related to the lowest eigenvalue of 
an appropriate equation, which is in turn, more or less, 
explicitly or implicitly, related to the important concept 
of energy quantization 0, 0, El flil h~f| . This traditional 
framework of determining upper critical fields by energy- 
quantization related techniques or scenarios has existed 
for fifty years or so 0,00- However, in this Letter, a 
new theory will be reported and justified to determine 
upper critical fields without recourse to energy quantiza- 
tion. Some applications and implications of the theory 



proposed will also be addressed. 

Various superconductors are layered compounds al- 
though they may be oxides, metallic or organic. To de- 
scribe layered superconductors, a continuous Ginzburg- 
Landau (CGL) model was proposed fl6|. in which the 
GL coefficients and the superpair masses are assumed 
spatially dependent. This model was considered to 
approach the limiting cases of the anisotropic Ginzburg- 
Landau (AGL) theory (see Ref. and references 
therein) and the Lawrence-Doniach (LD) model [lij . 

Recently, a modified CGL model has been proposed 
flij. Eof . The unit cell in this model consists of alternat- 
ing superconducting and weakly superconducting layers. 
The z-axis is normal to the layers and its origin is at 
the midpoint of one of the weakly superconducting lay- 
ers. The center of the superconducting layer is located 
at D/2, where D is the size of the unit cell [IsL l2l| . 
Applying an external magnetic field B parallel to the 
layers, the linearized CGL equation may be written as 
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where the nucleation center (z c — hk x /2eB ) is set at 
D/2. The condensation coefficient a(T,z) and the effec- 
tive masses, M(z) and m(z), are assumed as [Tol Eft Elj] 



a(T, z) = [a + at cos(2nz/D)} (1 - T/T c ), (2a) 
' G + Gicos(27rz/L>), (2b) 
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Here a>o, cui, Go, Gi, go and g\ are model parameters 
and the determinations of their values can be found in 
Ref. H3- 

To determine the upper critical field, Eq. Q should be 
completed by appropriate boundary conditions. Here we 
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choose the open boundary conditions (OBCs), 
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Eqs. (JU and (SJ can be transformed into a system of the 
form 



[/* =0, 



(4) 



where ^ is a wave function representing the discrete so- 
lutions of Eq. and U is the corresponding coefficient 
matrix j^. For Eq. Q to have non-trivial solutions, the 
determinant of U should be zero, 



det\U\ = 0. 



(5) 



It can be verified |22| that the parameters of the magnetic 
field B appear only in the main diagonal of U so that we 
have 



det \P - B 2 I\ = 0, 



(0) 



where J is a unitary matrix and P is a sparse matrix inde- 
pendent of B. Thus, the maximum B, namely the upper 
critical field, can be deduced from the largest eigenvalue 
of the following eigen equation 



(7) 



where x is the eigenfunction of P and this function is an 
auxiliary field to Eq. ©. Having determined B C 2, one 
can obtain the corresponding order parameter by substi- 
tuting B C 2 back into Eq. (0J. 

It should be stressed that Eq. £p) can be written as 



Q%(z) = B 2 f{?. 



(8) 



where Q is an operator. Eq. (JHJ), together with the bound- 
ary conditions of Eq. J2Lcan be written as the following 
matrix eigen equation [22| , 



= B *lf, 



(9) 



where \& possesses the same meaning as in Eq. flfl. Q is 
the same as P in Eq. J7J). Hence, Eq. © is equivalent to 
Eq. @ and *S? plays the same role as x as eigenfunctions. 

For comparison of our determinations with the tra- 
ditional treatment for bulk superconductors El H El 
l24j . we model the MgB 2 superconductor J25| whose 
anisotropy is small [20I Eq . For convenience, we treat 
it to be isotropic and ignore the spatial dependences in 
Eq. (0). Hence, Go = go = 1/m = 0.5 and G\ = gi = 
(Xi = were chosen and the nucleation center set at 
z c = 0. Other parameters such as an for MgB2 are taken 
from Ref. [2jj. It was found that (see Fig.QJ the results 
obtained from Eq. © are strikingly in agreement with 
the traditional GL work [III III El], 



B c2 = $ /2< 2 = V2kB c , 



(10) 
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FIG. 1: Temperature dependences of the upper critical field 
of MgB2. The calculated data (•) are in excellent agreement 
with the traditional GL theory of Eq. HUH (□) and reasonably 
consistent with the experiments in Ref. |26j [data (O) from 
the magnetization measurements; (A) and (v) from the mea- 
surements of the resistance-temperature curve near the super- 
conducting and normal states, respectively]. For comparison, 
results from the BCS theory and the two fluid model (2^1 are 
also illustrated. 



where the flux quantum $0 = h/2e and the coherence 
length £ = h/ y/2m\a(T)\ = h/^/2\a \(l - T/T c )/g . B c 
is the thermal dynamic critical field and K is the ratio 
of the penetration depth A to the coherence length £. 
It should be emphasized that Eq. (|10f> is derived from 
energy quantization of the quantum harmonic oscillator 
equation Q, H |24|] and the microscopic theory [f| may 
also arrive at a similar expression. The agreement of 
our calculations with Eq. (|10[1 lends us strong credence 
to our theory of upper critical fields. Furthermore, we 
found that the calculated results are reasonably consis- 
tent with the experiments of Finnemore et al. |2rj| and 
with some theoretical results 22, 27], as shown in Fig.^ 
These agreement and consistencies show that MgB2 is a 
superconductor describable within the GL/BCS frame- 
work. Note, however, that our main interest here lies in 
the qualitative predictions of the theories and thus we 
have assumed that the theories arrive at the same value 
of B C 2 at zero temperature; otherwise, the actual val- 
ues are deviated from the zero-temperature one but the 
qualitative trends remained. 

In Fig. |2 we have plotted the spatial distribution of 
the order parameter at B c2l which was obtained from 
Eq. |© at t = (T/T c ) = 0.9 for bulk MgB 2 . It can be 
seen that the order parameter is in excellent agreement 
with the following function at n = (in a.u.) 
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where Q = 1/ 1 ^/2B C 2 and H n are Hermite polynomials. 
Eq. (|llf) represents the eigenstates of the harmonic oscil- 
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FIG. 2: Order parameters at the upper critical field, the sec- 
ond and third largest square of the magnetic field respectively 
correspond to the ground, first and second excited states of 
the quantum harmonic oscillator [Eq. HH l. 



lator equation with energy serving as the eigenvalues [28j . 
The function at n — is the state of the lowest energy 
eigenvalue. The remarkable agreements of our calcula- 
tions with Eqs. i|10|) and unambiguously justifies the 
current theory of upper critical fields. 

It is interesting to find that (see Fig. |5J) the order pa- 
rameters corresponding to the second and third largest 
B 2 are nothing else but the first and second excited states 
of the harmonic oscillator equation, as given by Eq. l|Tl"|) 
at n — 1,2, respectively. Note that in the expression of 
C = 1/V2i?c2, B C 2 should now be replaced by B 2 and B 3 , 
which are the square roots of the second and third largest 
B 2 , respectively. Here, a few points are worthy of notice, 
(i) What does B associated with the lowest eigenvalue of 
B 2 mean? Which state are superconductors in subject 
to such B? The respective possible answers might be 
the lower critical field and the Meissner state? (ii) Are 
there any physical consequences in Eq. © [or Eq. 0] by 
changing 2e to e in Eq. Q? Are there any other phys- 
ical significances in Eqs. (JZJ) and © themselves except 
being used to obtain upper critical fields? (iii) The sign 
of ao in Eq. (|2a|) may need attention: by setting positive 
condensation energy «o (with a\ = 0), it was found that 
the calculated upper critical field is always zero. Hence, 
negative condensation energy is a necessity for a super- 
conducting state, as expected. 

Now some applications and implications of our theory 
will be addressed. The upward (positive) curvature of 
the B C 2-T curve in layered superconducting s yste ms is a 
subject of long-term interest H H E GJ, Hi M H |M 
El IHl IH Hill HI. In Fig.CHa), we have plotted the 
calculated parallel upper critical field as a function of the 
reduced temperature for the highly anisotropic cuprate 
superconductor Bi2Si2CaCu20s (Bi2212). The interval 
of the reduced temperature At for the plot is 0.0002. It 
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FIG. 3: Upper critical field and order parameter of Bi2212 
near T c . (a) Upper critical field as a function of reduced 
temperature t = T/T c , showing an upward curvature at t* ~ 
0.998. (b) Spatial distributions of the order parameters near 
T c . For t < t* , the order parameters are mainly confined 
within one unit cell whereas for t > t* , the order parameters 
spread over several unit cells. This crossover of the order 
parameter manifests itself as the upward curvature observed 
in (a). 



is clear that at t* ~ 0.998, an upward curvature appears. 
In the extreme vicinity of T c , the B C 2-T plot is linear 
(which is consistent with the GL and AGL theories ). A 
linear behavior also holds at low temperatures |2lj, |22| . 
Hence, a transition exists linking the two kinds of lin- 
ear behaviors. In Fig. EJ^b) , the spatial distributions of 
the corresponding order parameters are presented, where 
At = 0.0004 was chosen for clarity. It is clear that the or- 
der parameters below t* ~ 0.998 mainly reside within one 
unit cell [i.e., two dimensional (2D) behavior] whereas 
those above t* spread into the neighboring cells. With 
increasing temperature, more layers are penetrated [i.e., 
three dimensional (3D) behavior]. Thus, t* is a crossover 
temperature for the order parameter to transit between 
the 2D and 3D behaviors. Clearly, the upward curvature 
[Fig. Eta)] is originated from the crossover of the order 
parameter [Fig. E^b)], which is in turn due to that the 
layering structures [cf. Eq. @] come into effects. Note 
that the value of t* ~ 0.998 obtained here for Bi2212 is 
consistent with Ref. [13] , m which the most likely value 
was estimated in the range 0.997-0.998. 

It is worth mentioning that by considering the small 
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anisotropy in MgB2 (Go ~ 1, go ~ 1 and Gi = g± = ax = 
0), our calculations are consistent with the AGL theory. 
Moreover, CGL simulations have also been performed for 
MgB 2 but no obvious upward curvature was observed. In 
fact , th e issue of the upward feature is quite complicated 
HH 13 HH • A recent experimental study shows that 
the upward curvature may also exist in MgB2- If the 
same reason for the curvature, i.e., the crossover of the or- 
der parameter, is also applicable to MgB 2 , one can expect 
that the influence of layered structure on MgB2 should be 
less prominent than that on the highly anisotropic super- 
conductor Bi2212. Note that besides the upward feature 
obtained, we have also achieved a square-root field tem- 
perature dependence |2fl \2l\ within our framework. 



The theory reported here is generic. It can be used 
to study various properties of up per critical fields not 
only in different superconductors |20|, different bound- 
ary conditions [2ll 122 ] and different dimensions ; but 
also possibly in surface [22i and d-wave plsuperconduc- 
tivity and in the Gor'kov gap equation 6|. Moreover, 
the scheme in our theory may be used to determine the 
transition temperature of a superconductor in a magnetic 
field since the relationship between the upper critical field 
and temperature is equivalent to that between the critical 
temperature and the applied field. 



Our theory is an alternative to the traditional 
quantum-mechanical determinations of upper critical 
fields. The latter techniques are somehow related to the 
lowest eigenvalue (energy quantization) of the quantum 
harmonic oscillator equation. However, we have shown 
in this Letter that, instead of resorting to traditional 
method of energy quantization, the quantum harmonic 
oscillator equation can be treated just in a "classical" 
manner (i.e., directly utilize the largest eigenvalue of the 
square of the magnetic field) to determine upper criti- 
cal fields. It is known that some fundamental enigmas 
of quantum theory remain unresolved [40L 1-1 1 1 ] and tested 
[iij ]. Our success tempts us to ponder more about the 
techniques of quantum mechanics and to ask questions 
such as what is the origin of the wave function (which is 
the central puzzle of quantum mechanics [40()? It is clear 
that the concept of energy quantization can be bypassed 
in the present theory of upper critical fields. Moreover, 
the \& function in Eq. Q is equivalent to the x function 
in Eq. (JJJ) while % is an auxiliary field to Eq. JHJ). Hence VE* 
is an auxiliary field? Whether quantum mechanics needs 



interpretations or not 
mathematical beauty 
cess ji^. 



, the theory itself is of profound 
and has achieved brilliant suc- 
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